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A VARIATIONAL FORMULATION OF RELATIVISTIC HYDRODYNAMICS 

Nikodem J. Poplawski 

Department of Physics, Indiana University, Swain Hall West, 
727 East Third Street, Bloommgton, IN 47405, US45 

We combine Taub's and Ray's variational approaches to relativistic hydrodynamics of perfect 
fluids into another simple formulation. 
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The action for a relativistic perfect fluid in a gravitational field is given by the integral [l[ 

S=-J d±x^(-±R-pc 2 -pe-~ t i i (g^u»-l)+^\ (1) 



where R is the Ricci scalar, g is the determinant of the metric tensor g^ u , p is the rest density of the fluid (the number 



5S=- c J d^xy/^g^—^R^ - ^Rg^&gT ~ ftv"( fe ")^ -(c 2 + e + ^j5p 



qt< of particles per unit volume in the mean frame of reference of these particles), e(p, a) is the internal energy density 
per unit mass of the fluid, u M = is the four- velocity vector tangent to fluid-particle paths parametrized by A, a is 
the entropy density of the fluid, and p\ and P2 are Lagrange multipliers that ensure the normalization of and the 
l—— conservation of entropy for each particle, respectively. The variation of the action ([1]) is thus 

(N ■ 
> 

m : 

00 . 1 „ 1 \ 

(N ■ -pTScr + -{pc 2 + p^g Sg^ + -piu^ujg"" + ft,u" {8a) )(tl , (2) 

Csj; 2 2 / 

where we use the first law of thermodynamics, de = Tda + \dp (T - temperature, p - pressure), and the identity 
Su^ — u v {5x^) w . The semicolon denotes covariant differentiation with respect to the Christoffel symbols and the 
\ comma denotes partial differentiation. Integrating by parts in Eq. (|2|) gives 

■ SS= - c J d\^(^-^(R^ - ^Rg^Sg^ + {p lUl ,u v ),Jx^ - (c 2 + e + ^)5p 

-pTScr + ^(pc 2 + pe)g^Sg^ + ^p lUfl u v Sg^ - (^u^Sa^j . (3) 

In the integral ([TJ) we do not write explicitly, as in [l|, a Lagrange- multiplier term that corresponds to the conser- 
vation of particle-number density: 

(pa")* = 0. (4) 

Instead, we impose this condition through a relation between dp, Sg^ 11 and Sx^ in Eq. © 0. The equation of 
continuity (g]) gives 6 J d 4 x(y rl gpu ,1 )_ fl = or 5{^/~ r gpu tJ -) = 0, from which we obtain 

Sp=\pg lt »SgV-pupu''(6x*). v . (5) 
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Substituting Eq. © to © gives 



6S= - c J d^x^—g^-^R^ - \Rg^)5g^ - (pT + fau^S* 
+ ((/ii - pc 2 -pe- p^u"^ Jx^ - ^pg^Sg^ + ^lu^ujg^ . 



(6) 



Due to Hamilton's principle of least action, the action ^ vanishes for arbitrary variations 5a, Sx^ and 5g^ v . As a 
result we obtain, with the use of Eq. (TJ| , 

T , (7) 

\ p J ,u 

u v (-u^\ =0, (8) 

\ P / \v 

R^v - TjRg^ = k^u^Uv - pg^), (9) 

where 

A = m- pc 2 - pe-p. (10) 

The contracted Bianchi identity applied to the Einstein equations © gives, together with Eq. (U|), 

(pc 2 +pe+p + A)u tf u".„ = p,„(g^ - u»u v ). (11) 

If A ^ then Eqs. ((SJ) and (JTTJ) are in general not compatible because they constitute 8 constraints for only 5 
independent quantities: 3 components of w^ 1 , A, and e or p (which are related by the equation of state of the fluid). 
The compatibility of Eqs. |8]) and (fTTj) requires A = 0; Eq. © is thus satisfied identically and Eq. (fTTj) gives 4 
constraints for 4 independent quantities: 3 components of and e (or p). 1 Substituting the compatibility condition 
A = to Eq. (fTTj) gives relativistic Euler's equation of motion [HQ: 

(pc 2 +pe+ p)u v u^.^ = p !V (g^ - u»u v ), (12) 

which, in the limit c — * oo, reproduces the continuity and Euler's equation in nonrelativistic fluid mechanics 0]. 

Relativistic hydrodynamics of perfect fluids can be formulated from a variational principle with relativistic Clebsch's 
velocity potentials as dynamical variables 0] , but the presented formulation, like the methods of Taub and Ray, requires 
a smaller number of independent field components. Taub's formulation uses the Hclmholtz free energy density in the 
action and imposes the conservation of particle-number density through a relation between the variations of Sp, 8g^ v 
and fe M 0, [![, while Ray's formulation uses the internal energy density in the action and imposes the conservation 
of particle- number density through a Lagrange- multiplier term [l], Q . There also exists a Hamiltonian formulation of 
relativistic perfect fluids based on scalar fields, which is free of any constraints Q- 

The difference between the presented formulation and the variational approaches of Taub and Ray is that we use in 
the action (fTj) the internal energy density e, as in [l[, and use the variational relation ([5]) to impose the conservation of 
particle-number density, as in [2j. Because of this relation, we do not need to introduce Lin's constraint [l|,S|, without 
which isentropic flows described by the action for a perfect fluid with a particle-conservation Lagrange multiplier are 
irrotational *9|. We also, as in Ray's formulation, impose the conservation of entropy through a Lagrange- multiplier 
term in the action. In Taub's formulation, this conservation is derived from the variation of the quantity a related to 
T via We vary the particle paths, Sx^ (Lagrangian picture) Q, and not the four- velocity, 5u^ (Eulerian 

picture) [H liofl. although both approaches are equivalent in the presented formulation because we do not regard Sx^ 
as the coordinate variation and thus g^ v does not change under Sx^ . The metric tensor changes in Q under Sx^ and 
the corresponding field equation is the conservation of the dynamical energy-momentum tensor [ll| . 



For the equation of state of a dust, p = 0, Eqs. J5J and JTTJl are compatible. In this case, Eq. JTTJ reduces to the geodesic equation 
u v u fl . l/ = 0, which, with Eq. J5J, gives u u ^— J = J = 0- Consequently, the presented variational principle has a limitation for 

this case because A is not determined by the equation of motion; it is only constrained such that the ratio — is conserved for each fluid 
particle (including A = as a solution). 
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Our variational formulation of relativistic hydrodynamics of perfect fluids can be generalized to the Einstein-Cartan 
spacetime 12|, and to the general metric-affine spacetime |l3|. In the presence of a general affine connection, the 
continuity relation for a fluid takes its general-relativistic form (with the Christoffel symbols in the covariant derivative) 
O uM > otherwise the field equation resulting from the variation of the connection would impose unphysical constraints 



on the four- velocity b| . This field equation brings the connection to the Levi-Civita form combined with the torsion 



trace which does not appear in the symmetric part of the Ricci tensor so the Einstein equations remain unchanged. 
Therefore the presented formulation is valid also in the spacetime with a general affine connection, that is, with torsion 
and nonmetricity [15j. It should not be difficult to generalize this formulation to hyperfluids, which are continuous 
media that carry momentum and hypermomentum [161 ]. 
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